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ABSTRACT 

We present the 0.5 - 50 keV Suzaku broad-band X-ray spectral study of the Compton- 
thick AGN in NGC 5135. The Suzaku observation provides the first detection of NGC 
5135 above 10 keV that allowed us, for the first time, to estimate the absorbing column 
density, the intrinsic X-ray luminosity, the strength of the reflection component and 
the viewing angle of the torus for this AGN. The 0.5 - 10 keV spectrum of NGC 5135 
is characterized by the standard components for a Compton-thick source: a scattered 
continuum, a prominent Fe Ka emission line (EW ~ 2.1 keV) and a soft excess. At 
higher energies (E > 10 keV) the intrinsic AGN continuum shows up, implying an 
absorbing column density of the order of ^ 2.5 x 10^^ cm~^ and the intrinsic 2.0 - 10 
keV X-ray luminosity of ~ 1.8 x 10"*^ erg s~^. Assuming a toroidal geometry of the 
reprocessing material we show that an edge-on view of the obscuring torus is preferred 
in this source. 
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1 INTRODUCTION 

Seyfert galaxies host active galactic nuclei (AGN) pow- 
ered by accretion on to supermassive black hole. The lu- 
minous accreting engine is supposedly surroun ded by ob- 
scuri ng material having a torus-like geometry (|Antonuccil 
Il993l ). When the obscuring torus intercepts the observer's 
line-of-sight {i.e., type 2 AGN), the absorbing column den- 
sity (Nh) for X-ray photons originating from A GN is typi- 
cally higher than 10^^ cm"^ (|Cappi et al.ll2006l ). In mildly 
Compton-thick AGN (Nh ~ a few times of 10^"* cm"^), 
the primary emission is strongly suppressed below 10 keV 
and is seen only above 10 keV, while in heavily Compton- 
thick AGN (Nh ^ 10^^ cm~^) the primary emission is 
strongly depressed even at energies above 1 keV due to 
Compton down-scattering (|Matt et al.l |2000| ) . In Chandra, 
XMM-Newton and earlier observations sensitive only up to 
10 keV, Compton-thick AGNs have been identified only by 
using diagnostic properties such as high equivale n t wid th 
(EW) of Fe Ka line (> 1.0 keV, iBassani et all l|l999l) ). 
low flux ratio of har d X-ray (2.0 - 10.0 keV) to [GUI] 
A5007A line emission (jMaiolino et al.lll998l ) and there are 
no measurements of the true value of the absorbing col- 
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umn density. There is evidence that obscured AGNs are 
much more numerous than unobscured AGNs, both in 
the local universe a nd at intermed iate to high redshifts 
(IRisaliti et al. 1 1 19991 : iHasingeil I 20081 ). In recent years, ob- 
servations from INTEGRAL, Swift, and Suzaku have in- 
creased the number of C ompton-thick sources detected in 
hard X-ray jDella CecTe t al. 2008; Bcckmarm ct al. 20(51; 



iTueller et al 



20081 : ISevergnini et all l201lf ). However, the 



number of Compton-thick sources whose broad-band spec- 
tra are analyzed in detail is still limited. In this paper 
we present the Suzaku broad-band X-ray spectral proper- 
ties of NGC 5135, which is one of the brightest known 
Compton-thick Seyfert type 2 galaxies, with no previous X- 
ray observations above 10 keV. NGC 5135 has not been 
detected above 10 keV by any of the hard X-ray detec- 
tors, e.g., BeppoSAX-PDS, Swtft-BAT, INTEGRAL-IBIS, 
other than Suzaku HXD-PIN. NGC 5135 is a relatively 
nearby (redshift z ~ 0.014) galaxy and optically classified 
as a Seyfert t ype 2 on the basis of emission line ratios 
l|Phillips et al.1 ,1983). IR and UV studies have shown that 
NGC 5135 also contains a powerful compact nuclear star - 
burst (|Bedregal et al ] |2009l : iGonzalez Delgado et al.lll998l ). 
The higher spatial resolution of Chandra (~ 0.5") enables 
isolation of the AGN and starburst emissions, and has shown 
that the AGN in NGC 5135 is completely obscured by a col- 
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umn density Nh ^ 10^* cm"^ (|Levenson et al.ll2004 ). The 
Compton-thick obscuration around A GN in NGC 5135 i s 
also inferred from AS CA observations ijTurner et al.l[l997l ). 
iFukazawa et al.l (|201ll ) used the Suzaku data and report the 
presence of high absorbing column density (Nh) ~ 2.0 x 
10'^* cm~^ and a strong (EW ~ 1.6 keV) Fe Kg emi s sion l ine 
in NGC 5135. The primary aim of iFukazawa et al.l l|201ll ') is 
to study the properties of the Fe K line features and its 
correlation with the absorbing column density and AGN lu- 
minosity for a sample of Seyfert galaxies. In this paper, we 
discuss the detailed broad-band X-ray spectral analysis of 
NGC 5135 using the same Suzaku data, with the aim to 
characterize the reprocessing material around the AGN. 



by any other source. In order to do this, we searched the 
HEASARC/NED database for the potential hard X-ray (E 
> 10 keV) emitting sources in a circular region of the radius 
of 17' around the NGC 5135. The Swift, INTEGRAL and 
BeppoSAX catalogs do not report any hard X-ray source 
detection in this area. The Chandra, ASCA master catalogs 
report only NGC 5135 as the potential hard X-ray emit- 
ting source. During the Suzaku observations, NGC 5135 was 
placed at the HXD nominal pointing, therefore, in the spec- 
tral analysis, we used a cross-calibration constant of 1.18 
between the HXD and XIS spectra, as suggested by the 
Suzaku-YiXD calibration team. 



2 OBSERVATIONS AND DATA REDUCTION 

Suzaku observed NGC 5135 on 2007 July 03 (observation 
ID 702005010) with an exposure time of ~ 52.5 ks. Suzaku 
llMitsuda et all 120071') carries four X-ray telescopes (XRTs: 
ISerlemitsos et al.l 1,200?)) with X-ray CCD cameras (XIS) at 
their focal-planes. The XISs are sensitive to 0.2 - 12.0 keV 
energy band with 18' x 18' field-of-view. Among the four XIS 
CCDs, three (XIS 0, 2, and 3) are fron t-illuminated (FI) ari d 
one (XISl) is back-illuminated (BI) (|Kovama et al.ll2007h . 
There are no observations with XIS2, due to a malfunction 
in November 2006. Suzaku also has a non- imaging hard X- 
ray detector fHXD: pTakahashi et al.1 (120071 ) 1 which has two 
types of detectors, the PIN and the GSO, which are sensi- 
tive in 10 - 700 keV energy band. 

NGC 5135 was detected by XIS 0, 1, 3 and HXD-PIN. We 
obtained the XIS and HXD event files and reduced these by 
following the standard procedure described in the Suzaku re- 
duction guidon, and using the most recent calibration files. 
For the low-energy instruments (XISs), the source spectra 
were extracted from a circular region of 2.4' radius centered 
on the source. The background spectra were extracted from 
two circular regions of 2.2' radius offset from the source and 
calibration sources. The XIS response (rmfs) and ancillary 
response (arfs) files were produced with the ftools tasks xis- 
rmfgen and xissimarfgen, respectively, and the latest calibra- 
tion files were used. The spectra of the two front-illuminated 
CCDs (XISO and XIS3) were merged. The net XIS source 
spectra were binned to have minimum signal-to-noise ratio 
(S/N) of 4 in each energy bin and statistics has been 
used. For the HXD-PIN, we used the rev2 data that in- 
clude all 4 clu ster units and the best available background 
IFukazawa et al. 2009), w hich accounts for the instrumental 
background (NXB. lTakahashi et al.l (|2007l ')) with systematic 
uncertainties of ~ 1.3% (at 1 a). We then simulated a spec- 
trum for the cosmic X-ray background counts and added it 
to the instrumental one. Using this background NGC 5135 
is detected in the 13 - 50 keV band at ^ 7% above the back- 
ground with a net count rate of (2.7±0.2) x 10~^ cts s~^ (a 
total of ~ 1400 net source counts have been collected), cor- 
responding to a signal-to-noise ratio (S/N) ~ 9.2. Since the 
HXD-PIN is a non-imaging detector and has large field-of- 
view {r^ 0.56 deg X 0.56 deg), we checked that the HXD-PIN 
detection of NGC 5135 is not caused and/or contaminated 

^ http: / /heasarc. gsfc.nasa.gov/docs / suzaku / analysis / abc / 



3 SUZAKU X-RAY SPECTRAL FIT 

The 2.0 - 10 keV Suzaku data can be fitted with a flat power 
law (r ~ 1.4) and a Gaussian for Fe Ka line of high EW (~ 
3.0 keV), which is broadly si milar to the fit report ed by pre- 
vious Chandra observation (|Levenson et al. 1 120041 '). The flat 
power law and high EW of Fe Ka line are indicative of the 
heavily obscured AGN. In order to characterize the X-ray 
emission from heavily absorbed AGN, we attempt to flt the 

0. 5 - 50 keV Suzaku data using physically motivated models. 
In order to avoid the complexity of the soft component (E 
< 2.0 keV) we first attempted to flt the 2.0 - 50 keV spec- 
trum with a basic model of obscured AGN i.e., an absorbed 
power law representing the AGN X-ray emission transmitted 
through a cold absorber and an unabsorbed power law for 
the scattered component. This simple model gives x^/d.o.f. 
~ 182/31 and leaves large residuals around 6.4 keV with an 
emission line-like shape. Adding a Gaussian proflle for the 
Fe Kq at 6.4 keV improves the fit very significantly and gives 
the fit statistics x^/d.o.f. ~ 39/29, with Nh ~ 2.5 x 10^* 
cm~^, F ~ 1.8 for the absorbed component, a steep power 
law (r ~ 2.7) accounting for the emission below 10 keV and 
the equivalent width (EW) of Fe Ka line ~ 2.1 keV. The 
Fe Ka line is fitted with a narrow Gaussian {a fixed to 10 
eV) and any increase in the line width worsen the fit. There 
are no signatures for the Fe K/3 or the ionized component 
of Fe Ka line emission in the residuals and the addition of 
any Gaussian component for such lines worsen the fit statis- 
tics. If the photon index (F) of the unabsorbed power law 
component is fixed to the value equal to the photon index 
of the absorbed power law component, the fit statistics does 
not improve and this fit results F ~ 2.1 with little increase 
in the column density ( Nh ~ 2.6 x lO'^* cm~'^), while the 
other parameters remain nearly unchanged. Considering the 
high Nh of the absorbed component and the high EW of 
the Fe Ka line as the indicative of t he presence of reflection 
comp onent, we added the PEXRAV IjMagdziarz fc Zdziarskil 
I1995I ) model component which represents the emission re- 
flected from a neutral medium. The addition of reflection 
component gives x^/d.o.f. ~ 33/28 (model 'Ml' in Table 1) 

1. e., improvement of A^^ ~ 6 for 1 d.o.f. at the signiflcance 
level of 90%. The photon index of the reflected component 
(PEXRAV) was fixed equal to the transmitted component and 
the reflection scaling factor was fixed to -1, while all other 
parameters of the PEXRAV model were kept to their default 
values. We note that an equally good fit is obtained if we 
fix the photon-index of the intrinsic power law to 1 .9, the 
canonical value for Seyferts l|Nandra fc Pouiid^l 19941 ). How- 
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ever, as expected, a higher value of the absorbing column 
density (Nh ~ 2.9 x 10^'' cm~^) is required to account for 
the steeper power law. In order to assess the importance 
of the reflection component, we estimated the 2.0 - 10 keV 
X-ray fluxes associated with the unabsorbed power law and 
the PEXRAV (reflection) model components using the model 
'Ml'. The 2.0 - 10 keV X-ray fluxes associated with the 
unabsorbed power law and the PEXRAV (reflection) model 
components are ~ 1.27 x 10"^'^ erg cm^^ s~^ and ~ 1.49 
X 10~^^ erg cm"'^ s~^, respectively. These flux values corre- 
spond to ^ 28% and ~ 34% of the total observed flux in 2.0 
- 10 keV band, respectively, and suggest that both the unab- 
sorbed power law and the reflection components contribute 
substantially in this band. 

Since Compton scattering is signifl cant for large col- 
umn densities and high energy photons (! Yaaooblll997l ). we 
applied the scattering component using CABS model to the 
primary transmitted component (model 'Ml^' in Table 1). 
However, the CABS model does not account for the photons 
scattered into the line-of-sight from other directions. To ac- 
count for this effect, we fitted the spectrum with a model 
('Ml*' in Table 1) in which the PLCAB model accounts for 
Compton scattering assuming a uniform, spherically dis- 
tributed regrocessing material around the X-ray emitting 
source (|Yaaoobl Il997l ). This model also gives fit statistics 
(x'^/d.o.f. ~ 32.5/27) and parameters similar to the models 
'Ml' and 'Ml*' (c/., Table 1). 

We also checked the spectral fitting with a model ('M2' in 
Table 1) in which the reflection component is accounted by 
the REFLIONX model (Ross & Fabian 2005). The REFLIONX 
model characterizes the reflected emission from an optically 
thick disc illuminated by radiation with a power law spec- 
trum and produces the florescent emission lines Eis well as 
the continuum emission. The best-fit which could fit the Fe 
Ka emission line gives unusually high Fe abundance (~ 10 
times Solar), low ionization parameter (^ ~ Ollyg^ erg cm 
s"^) and photon index (F) of ~ 1.6 with xV<i-o.f. ~ 37.6/28. 
The required high abundance and the low ionization param- 
eter suggest that the Fe line emission is originating mainly 
from neutral or mildly ionized material. 
Since the PLCAB model is valid up to energies between 10 and 
18.5 keV and for column densities up to ~ 5 x 10^** cm"'^ 
l Yaaooblll997h. we attem pted to use the MYTORUS model 
i Murphv fc Yaqoob|[200 9') that accounts for the photoelec- 
tric absorption, Compton down-scattering as well as Fe Ka, 
K/3 line emission and is valid for the 0.5 - 500 keV en- 
ergy band and for column density up to 10^^ cm~^. This 
model considers an azimuthally symmetric toroidal geom- 
etry around AGN for the X-ray reprocessor and produces 
the reprocessed continuum and Fe Ka, Fe K/3 lines self- 
consistently. The opening angle of the torus is fixed to 60° 
and the column density as well as the viewing angle are free 
parameters. The best-fit using the MYTORUS model ('M5') 
renders Nh ~ 2.7 x 10^* cm~^, F ~ 1.96, and the viewing 
angle of the torus ~ 90°, with xVd-O-f- ~ 45/27. The best- 
fit is obtained by fixing the viewing angle to 90° . If we allow 
it to vary, the fit statistics does not improve {i.e., similar 
value with decrease of 1 in d.o.f.) and gives viewing angle ^ 
88° at 90% confidence level. We note that while using the 
MYTORUS model ('M5'), an additional unabsorbed power law 
of photon index ~ 2.0 is required to flt the data below 10 
keV. 
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Figure 2. Confidence contours for the photon index versus ab- 
sorbing column density obtained from model 'M3'. Contours at 
1(7, 2a, and 3cr level are shown in Black, Red and Green, respec- 
tively. 

We flnally fitted the whole 0.5 - 50 keV spectrum with the 
above models and noted that addition of a thermal plasma 
model (MEKAL in XSPEC) of temperature kT ~ 0.70 keV, is 
required to account for the emission below 2.0 keV (mod- 
els 'M3', 'M3I"', 'M3*' and 'M4' in Table 1). While fitting 
the 0.5 - 50 keV spectrum, we have excluded the XISs data 
points lying in 1.6 - 1.9 keV energy range because they are 
likely to be affected by systematic calibration uncertainties 
around the instrumental silicon K-edge. We find that the 
best-fit parameters obtained are similar in both the 2.0 - 50 
keV and 0.5 - 50 keV spectral fittings. This implies that the 
addition or removal of the soft component (E < 2.0 keV) 
does not affect the interpretation of the hard components. 
Table 1 lists the models and their resulting best-fit param- 
eters. All the quoted errors are at 90% confidence level for 
one interesting parameter. The best-fitted spectrum using 
model 'MS' and the unfolded model are shown in Figure 1. 
Figure 2 shows the confidence contours for the photon index 
versus absorbing column density using model 'M3'. 



4 DISCUSSION 

4.1 Soft X-ray emission 

In our Suzaku X-ray spectral fltting, the soft X-ray emission 
below 2.0 keV is accounted by a thermal plasma model 
of temperature kT ~ 0.70 keV and a steep unabsorbed 
power law. The observed 0.5 - 2.0 keV luminosity is ~ 1.66 
X 10*^ erg s~^ (c/., Table 2), consistent with previous 
ASCA and Chandra observations within the measurement 
uncertainties. The ratio of soft X-ray to far infrared 
emission (Lo. 5-2.0 kov/Lpm) ~ 3.0 x 10"'' for NGC 5135 
is consistent w ith the typica l ratio observed in starburst 
galaxies {e.g., *Ra nalli et al.l ()2003l )l. suggesting that the 
soft X-ray emission is dominated by starburst emission. 
Indeed, Chandra observation of high spatial resolution 
shows that the 0.5 - 2.0 keV soft X-ray emission from 
the central region is dominated by the circumnuclear 
starburst. The soft X-ray spectral components i.e., thermal 
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Table 1. The best-fit spectral parameters 



Energy 


Model 


Soft Component 


Hard Component 








range 




kT 

(kcV) 


Tsoft 


Nh 

(x 10^4 cm^2-) 




R 


Elinc 

(kcV) 


EWpc 
(kcV) 




2.0 - 50.0 


M5 




^•"^-0.08 


9 71+0.81 
^- '^-0.75 


1 96+0.26 
^■^"-0.08 








44.9/27 


2.0 - 50.0 


Ml 




^•'"-0.25 


9 r:ri+1.28 


1 ^o+io'^ 

i.OU_o 7j 


^■^'^-0.27 




9 01+0.44 
^•^^-0.27 


32.8/28 


2.0 - 50.0 


Mlt 




'7C+O.22 
^•'"-0.18 


9 CI +1.06 
^•^^-1.08 


1 

i-Oi_0.66 




6.38«:«1 


9 14+0. 36 
^•^^-0.38 


32.6/28 


2.0 - 50.0 


Ml* 




^■»»_0.22 


9 00+1. 03 
^•■J'->_0.91 


i.68ir7? 






9 og+0.42 
^•^^^-0.38 


32.5/27 


2.0 - 50.0 


M2 




r, -17+0.78 
' -0.52 


9 04+0.56 
^•■J^-0.36 










37.6/28 


0.5 - 50.0 


M3 




r, 70+0.19 
^- '•^-0.17 


9 ^Q+1.04 


1 ^t; + 1.09 


^^■4^-0.21 




9 -,(,+0.46 
^•^^'-'-0.32 


89.0/70 


0.5 - 50.0 


MS"!" 


n 70+0.07 

'J- ''J-0.03 


9 09+0. 18 
^■°^-0.18 


9 74+1.76 
^•'^-1.05 


1 50+102 

l.O»_0 gg 




6.38tO;«l 


9 07+0.44 
^•'^'-0.38 


91.4/70 


0.5 - 50.0 


MS* 


fi7+0-08 
-0.03 


9 7-1 +0.19 
^- ' -"--0.18 


9 on+1.00 


1 fiO+''-^'^ 
i.OU_o 79 




6.38tO:«l 


9 r|Q+0.40 


89.0/70 


0.5 - 50.0 


M4 




9 01 +0.12 
^■•^-'--0.07 


9 9n+0.36 
^•^'^-0.27 


l-00_o 12 








107.7/72 



Notes:- Ml: const(pl-|-wabs*pl-|-pexrav-|-line), Ml^: const(pl+wabs*cabs*pl+pexrav+line), 
MlJ: const(pl+plcabs+pexrav+line), M2: const(pl+wabs*pl+refiionx.mod), 

M3: const(mckal+pl+wabs*pl+pexrav-|-line), M3^: const(mekal-(-pl-|-wabs*cabs*pl-|-pexrav+line), 
M3*: const(mekal-Hpl-|-pl cabs+pexrav+line), M4: const(mekal+pl+wabs*pl+reflionx.mod), M5: MYTORUS 
jMurphv fc Yaooobllioogh . 
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Figure 1. Suzaku 0.5 - 50 keV spectral fit using model 'M3'. The left panel shows the spectral fit and residuals (in the bottom panel) 
and the right panel shows the unfolded model 'M3'. The XISO-l-3, XISl and HXD-PIN data points and spectral components are shown 
in Black, Green and Red, respectively. In unfolded model, the additive components are shown by dotted curves and cumulative model is 
shown by solid curves. 



component of temperature kT ~ 0.70 keV and a steep 
power law (F ~ 2.7), in our Suzaku data are fairly consistent 
with the X-ray spectral properties of a spatially-resolved 
strong circumnuclcar starburst region identified in Chandra 
observation ( Levenson c t al. 2004). The steep unabsorbed 
power law is likely to have contribution from different 
individ ual sources, mainl y X-ray binaries in the starburst 
region jZezas et al.l |200"3 ). High resolution near-IR study 
of the central ~ 2.3 kpc region in NGC 5135 has revealed 
the star formation knots and a large ^ 600 pc [Si VI] 
1.96Atm line emissio n ionization cone centered on AGN 
l|Bedregal et al.ll2009l ). It has been argued that SNR shocks 
play a dominant role in ionizing the gas in central ~ 2.3 
kpc region and photoionized emission due to AGN as well 
as the recent star f ormation have rather localized and 
smaller contribution (|Bedregal et al.l |2009| ). We note the 
fact that in general, high resolution X-ray grating spectral 
studies show that the soft X-ray emission in Seyfert 2s 



is dominated by recombination lines implying photoion- 
izatio n as the primary mechan i sm llKinkhabwala et al.l 
l2002l : iGuainazzi fc Bianchil l2007l : iMarinucci et all 120111 ) 

However, the soft X-ray emission in Seyfert 2s with intense 
circumnuclear starburst, is likely to be dominated by coUi- 
sionally ionized opti cally thin plasma em ission associated 
with star formation (|Guainazzi et al.ll2009l ). Limited spatial 
and spectral resolution of the Suzaku XIS CCD does not 
allow us to quantify the AGN and starburst contribution to 
the 0.5 - 2.0 keV soft X-ray emission, which in turn, unable 
us to draw firm conclusions about the origin of the soft 
X-ray emission. X-ray observations of high spectral and 
spatial resolution are required to confirm the nature of the 
soft X-ray emission in NGC 5135. 



5 



Table 2. Suzaku X-ray fluxes and luminosities of NGC 5135 



Energy Model 
(keV) 



-^obs 

(erg cm~^ 



Lobs 

(erg s-i) 



0.5 - 2.0 


M3 


3.89 


X 


10" 


13 


1.66 


X 


10*1 


5.0 




MSt 


3.89 


X 


10- 


13 


1.66 


X 


10« 


65.9 




MS* 


3.89 


X 


10- 


13 


1.66 


X 


10« 


16.6 




M4 


3.88 


X 


10- 


13 


1.66 


X 


10« 


12.3 


2.0 - 10 


MS 


4.44 


X 


10- 


13 


1.89 


X 


10« 


11.3 




M3t 


4.40 


X 


10- 


13 


1.87 


X 


10« 


124.8 




MS* 


4.43 


X 


10- 


13 


1.89 


X 


10« 


27.0 




M4 


4.37 


X 


10- 


13 


1.86 


X 


lO^i 


20.5 




M5 


4.17 


X 


10- 


13 


1.78 


X 


1041 


135.5 


10 - 50 


MS 


1.58 


X 


10- 


11 


6.67 


X 


10*2 


1.2 




M3t 


1.50 


X 


10- 


11 


6.39 


X 


10-12 


8.4 




MS* 


1.50 


X 


10- 


11 


6.67 


X 


10*2 


1.7 




M4 


1.49 


X 


10- 


11 


6.35 


X 


1012 


1.3 




M5 


1.46 


X 


10- 


11 


6.17 


X 


1012 


4.5 


Fe Ka 


M3 


9.35 


X 


10- 


14 


3.98 


X 


lOlo 





Notes:- Fobs^ observed flux, Lobs^ observed luminosity, Fcor: 
absorption corrected flux. 



4.2 Hard X-ray emission 

Our Suzaku X-ray spectral fitting sliows tliat the 2.0 - 10 
IceV continuum emission has contribution from several spec- 
tral components e.g., steep unabsorbed power law, ther- 
mal plasma model, reflection component, absorbed power 
law (see, Figure 1). However, the relative contribution from 
each individual spectral component is energy dependent. 
For example, the relative contribution of the unabsorbed 
power law decreases towards higher energies, while it in- 
creases for the reflection component. The thermal plasma 
model and absorbed power law components contribute only 
to the soft and hard part of the 2.0 10 keV band, re- 
spectively. Using SFR - X-ray correlation (SF R ~ 2.0 x 
L2.0-10 kcv M0 yr~i) from lRanalli et ahl (|2003|) and assurn - 
ing SFR ~ 15 Mq yr"! in NGC 5135 (|Bedregal et a"l]|2009^ ■ 
we noted that the estimated contribution from the circum- 
nuclear starburst to the 2.0 - 10 keV X-ray luminosity is 
~ 7.5 X lO*" erg s'l, which ~ 40% of the total observed 



2.0 - 10 keV luminosity. This value is similar to the 2.0 
- 10 keV non-AGN X-ray luminosity reported in Chandra 
observation, which could spatially res olve the A GN and cir- 
cumnuclear starburst X-ray emission (|Levenson et al. 2004) . 
Above argument infers that nearly ~ 60% of the total ob- 
served 2.0 - 10 keV X-ray emission is attributed to AGN 
emission, which is accounted by the scattered, reflected and 
transmitted components. Our Suzaku (0.5 - 50 keV) broad- 
band X-ray spectral analysis of NGC 5135 shows the char- 
acteristics of starburst as well as obscured AGN and thus 
confirming its composite nature. Furthermore, the X-ray 
emission from circumnuclear starburst regions in NGC 5135 
can be characterized by both thermal plasm a model plus 
unabsorbed power law (|Levenson et al.l |2004| '1. It is noted 
that star-forming galaxies show characteristic power law X- 
ray emission that co uld be the integra ted spectrum of dis- 
tinct X-ray sources (|Zezas et all 120021 '). Moreover, an un- 
absorbed power law component representing the scattered 
emission of the primary AGN continuum to the line-of-sight 
is also commonly seen in obscured AGNs llBraito et al ]|2009l : 
IComastri eraLH2010l : ISevergnini et al.ll201ll 'l. With present 



Suzaku data we cannot spatially separate out the AGN and 
circumnuclear starburst emission and therefore it is possible 
that the unabsorbed power law component in our spectral 
fitting may have contribution from both the circumnuclear 
starburst as well as the scattered component of AGN emis- 
sion. It is likely that the quality of the present 2.0 - 50 keV 
Suzaku data do not allow us to accurately determine the rel- 
ative contributions from starburst, scattered, refiected and 
absorbed AGN emission spectral components in the 2.0 - 10 
keV energy band. 

Our Suzaku 2.0 - 50 keV broad-band spectral fitting 
clearly shows that the hard X-ray continuum is charac- 
terized by a heavily absorbed power law and a reflection 
component, con s istent with the previous results reported by 
iFukazawa et al.l (|201ll ). The high absorbing column density 
(Nh ~ 2.5 X lO^^* cm~^) and the high ratio of the intrin- 
sic to the observed 2.0 - 10 keV flux/luminosity (see. Table 
2) infer the presence of a heavily obscured, luminous AGN. 
The ratio of the intrinsic to the observed flux/luminosity in 
2.0 - 10 keV band, found for NGC 5135, is consistent with 
the typical values obs e rved in other Compton-thic k AGNs 
(|Levenson et al.l I2OO6I : iBrightman fc Nandral I2OIII ). Given 
the importance of Compton scattering in absorbers with col- 
umn densities of a few times of 10^* cm~^, the absorption 
corrected (intrinsic) flux/luminosity is model-dependent. 
The value of the intrinsic luminosity ( il2.o^_Ji±£v^ ^ 125) 

^ ^2.0 - 10 kcV.obs ^ 

obtained from the model 'M3^' which uses the CABS model 
for the scattering, can be considered as the upper limit. Since 
the CABS model represents the Compton scattering from the 
matter only along the line-of-sight and therefore offers a 
small covering fraction i.e., a single cloud along the line- 
of-sight which only removes photons from the beam owing 
to Compton scattering and does not consider the addition 
of photons scattered from other directions into the line-of- 
sight. Moreover, the ratio of the intrinsic to the observed 
flux/luminosity also depends on the photon index of the 
intrinsic power law, for example, this ratio obtained from 
M3t model is as high as ~ 300, if we fix the photon index to 
1.9. The value of the intrinsic luminosity obtained from the 
model 'M3"''' where the PLCABS model is used for scattering 
can be considered as a lower limit, since in this scenario, the 
contribution from photons initially not in the line-of-sight, 
which are then scattered to the line-of-sight, is highest, and 
roughly compensate the loss of photons which are initially 
along the line-of-sight, and that are scattered out. Thus, 
the two models represent two limits: a single cloud covering 
(CABS) gives an upper limit and full covering (PLCABS) gives 
a lower limit on the intrinsic luminosity, and the true value 
lies somewhere in between depending on the covering factor. 
In the MYTDRUS model covering factor corresponds to 0.5 and 
the ratio of the intrinsic to the observed flux/luminosity (~ 
135) indeed lies between the two limits, if the photon index 
of the intrinsic power law is assumed to be 1.9. 
High obscuration around the AGN is also inferred from the 
flux ratio of hard (2.0 - 10 ke V) X-ray to [OIII] A5007A line 
emission (|Bassani et al.lll999l 'l. i.e., for NGC 5135, we obtain 
log[F2.o - 10 keV, obs/F[oiii]] 'T ^ -1.1, usiug the reddeuiug- 
corrected [OIII] flux given in ISingh et al.l (|201lh . Further- 
more, the flux ratio of the intrinsic 2.0 - 10 keV hard X-ray 
to [OIII] A5007A line emission is (F2.0 - 10 kcV, int/F[oiii]) 
^ 9.2, consistent with the mean ratio of observed in Seyfert 
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1 galaxies l|Heckman et alfeoOSh . supporting the notion that 
the obscured nucleus in NGC 5135 is intrinsically similar to 
a Seyfert type 1 nuclei. 

Signature of Compton-thick matter around the AGN can 
be characterized by Compton reflection continuum with a 
broad hump peaking around 20 - 30 keV, which rapidly 
decreases at both low and high energies due to photoelec- 
tric absorption and Compton down-scattering, respectively. 
In our spectral fits the reflection strength measured as the 
ratio of the normalizations of the reflection component to 
the transmitted component (absorbed power law) is ~ 0.4 
(using model 'Ml'). We do not attempt to compare the 
strengths of the reflection components obtained by using 
the MYTORUS model and the PEXRAV model since these mod- 
els consider different geometries for the reprocessor. The 
relative strength of the reflection continuum very critically 
depends on geometry since it is affected by the angle of 
reflection integrated over the surface of the reprocessor. 
iMurphv fc YaaoobI (|2009l ) showed that the reflection spec- 
trum from a Compton-thick face-on torus that subtends the 
same solid angle as a Compton-thick face-on accretion disc 
at the X-ray source is much weaker and it is not meaning- 
ful to compare the reflection strength given by the MYTORUS 
model with one from the PEXRAV model. 

4.3 Geometry of the absorber 

It is important to note that the PEXRAV as well as REFLIONX 
models consider reflection from the surface of accretion disc 
and therefore, are insufficient to constrain the geometry of 
the toroidal-shaped obscuri ng material around the A GN. 
We used the MYTORUS model (|Murphv fc Yaqoob|[2009l ) that 
considers the geometry of reprocessing material as an az- 
imuthally symmetric obscuring torus. This model suggests 
the edge-on view (i.e., the viewing angle ~ 90°) of the ob- 
scuring torus. The absorbing column density given by the 
MYTORUS model corresponds to the equatorial column den- 
sity of the toroidal-shaped obscuring material. We note that 
the column densities obtained from the MYTORUS and other 
models are similar within the uncertainties (see, Table 1) 
and therefore, it again confirms that the obscuring torus is 
oriented nearly edge-on. We caution that in the current ver- 
sion of the MYTORUS model, the torus opening angle is fixed 
to 60° while it needs to be a varying parameter. The open- 
ing angle of 60° corresponds to covering factor of 0.5 and the 
solid angle subtended by the torus on to the X-ray emitting 
source to 2n. Hints on the covering fraction of the absorber, 
as seen from the X-ray source, can be inferred from our fits. 
It has been shown that variation in the geometry of ob- 
scuring torus, and not the iron abundance or intrinsic spec- 
tral shape, is required to produce EWs significantly larger 
than 1 keV and therefore, EW of Fe Kg line can be used t o 
place limits on the co vering factor (|Levenson et al.ll2002h . 
iLevenson et al.l l|2002l ') argued that the high EW of Fe Ka 
line (~ 2.1 keV) in NGC 5135 suggests a small opening angle 
(~ 20°) of the obscuring torus giving a covering fraction of 
~ 90%. Our Suzaku X-ray s pectral fitt i ng giv es Lpe ~ 3.98 x 
10*° ergs"^ and using Leve nson et al.l (|2006t ) empirical ratio 
of Fe Ka line luminosity to intrinsic 2.0 - 10 keV luminosity 
we obtain L2-W kov, int ~ 1-99 x 10*^ erg s~^, i.e., nearly 
100 times the observed 2.0 - 10.0 keV luminosity. This value 
of the 2.0 - 10 keV intrinsic luminosity lies between the val- 



ues from the MYTORUS model (covering fraction ~ 50%) and 
the MS* model that uses the PLCABS model (covering frac- 
tion ~ 100%). Furthermore, we note that the intrinsic 2.0 - 
10 keV X-ray luminosity obtained by using the mid-IR (12 
/im) - X-ray luminosity correlation (Lmir - L2.0-10 kcv) and 
mid-IR luminosity (Lmir) ~ 1-15 x lO'*'^ erg s~^ given in 
iGandhi et ai] Hooi), is ~ 7.6 x lO"^ erg s"\ that again lies 
between than the values from the model M5 and M3"'". The 
high resolution mid- IR core fluxes/luminosities presented in 
iGandhi et al] (|2009| ) are reported to be least contaminated 
by the circumnuclear starburst and primarily represent the 
emission from the torus. The above comparisons for the 2.0 
- 10 keV intrinsic luminosities suggest that the covering fac- 
tor lies between 50% to 100%. Finally, the low ratio between 
the soft power law flux and the intrinsic emission (< 0.01, 
Table 2) suggests a low scattering efficiency, which again, 
may be due to a high covering factor of the thick absorber, 
leaving only a small fraction of the solid angle free for the 
primary continuum to reach the warm scatterer. A similar 
low flux ratio (< 0.01) of the unabsorbed powerlaw to the 
intrinsic 2.0 - 10 keV emission is obtained with the model 
'M5' which considers an axi-symmetric toroidal geometry of 
the reprocessing material around the AGN. The actual scat- 
tered fraction may be even lower, given that part of the soft 
power law may be due to the starburst emission. 



5 CONCLUSIONS 

We analyzed the Suzaku XIS -I- HXD-PIN observation of the 
AGN in NGC 5135. We found that the 0.5 - 50 keV spec- 
trum can be reproduced by a model consisting of; (i) a soft 
component, characterized by a thermal plasma model (the 
MEKAL) of temperature kT ~ 0.70 keV plus an unabsorbed 
steep power law (F ~ 2.7), (ii) a hard component best fltted 
by an absorbed power law (Nh ~ 2.5 x 10^* cm~^, F ~ 1.6), 
a reflection continuum component and a prominent Fe Ka 
line with EW ~ 2.1 keV. Suzaku broad-band energy cover- 
age allowed us to accurately measure the absorbing column 
density, intrinsic AGN luminosity and the reflection com- 
ponent. We attempted various physically motivated models 
and all the models confirm that the AGN in NGC 5135 is 
obscured by Compton-thick material. The MYTORUS model 
which considers the geometry of the reprocessing material 
as a azimuthally symmetric torus suggests that the torus is 
viewed nearly edge on. Both the estimates on the reflection 
strength and the comparison between the emission line flux 
and the continuum flux suggest that the obscuring torus 
may cover 90% of the line-of-sight. Our Suzaku X-ray spec- 
tral study of NGC 5135 may represent a case study of some 
of the complex features in Compton-thick obscured AGNs. 
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